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We study the electronic transport in an infinite one-dimensional Hubbard chain, driven by a 
homogeneous electric field. The physical chain is coupled to fermion bath chains, in order to account 
for dissipation and to prevent the occurrence of Bloch Oscillations. The steady state current is 
computed in the frame of Keldysh Green’s functions in Cluster Perturbation Theory. The current 
characteristics are dominated by resonant-tunneling-like structures, which can be traced back to 
Wannier-Stark resonances due to anti-ferromagnetic correlations. The same current characteristic 
occurs in a non-interacting Wannier-Stark model with alternating on-site energies. Non-local effects 
of the self energy can be accounted for the observed physical behaviour. 

PACS numbers: 71.15.-m, 71.27+a, 71.10. -w, 73.23.-b 


I. INTRODUCTION 

One dimensional (Id) quantum systems exhibit fas- 
cinatiiM non-equilibrium effects such as Bloch oscil- 
lationd^ (BOs), Zener tunneling and Wannier-Stark 
resonances.l^lt^ The first theoretical studies of these ef¬ 
fects date back to the early days of condensed matter 
physics and were followed by long-lasting controversial 
discussions about their relevance in real materials.^^^ It 
was only within the last two decades that the theoret¬ 
ical predictions could be observed in condensed matter 
sytems and optical potentials in a series of ingenious 
experiments.!^*^. A detailed overview over the most 
important experimental verifications of the above men¬ 
tioned effects is provided in Ref. HO] and Ref. [HJ where 
also experimental set-ups for determining Wannier-Stark 
resonances with cold atoms in optical lattices are intro¬ 
duced. 

Bloch oscillations are mainly studied in Id systems 
which are and were frequently used as model systems to 
investigate the physical nature of BOs. As far as corre¬ 
lations are concerned the Id Hubbard model was mostly 
the model of choice. The clear evidence for the existence 
of BOs and Wannier-Stark resonances led to a series of 
further theoretical investigations of Id structures under 
the impact of a unifo rm force induced by an electric or 
magnetic field llZ120|22f(25| 

For the emergence of Wannier-Stark resonances in Id 
systems a periodic potential has to be applied. It was 
long thought that this periodic potential has to meet 
some specific conditions Recently, it was shown ana¬ 
lytically that any periodic potential suffices to cause these 
resonant quantum states.!^ 

The theoretical investigation of non-equilibrium prop¬ 
erties have become very popular in recent years stimu¬ 
lated by the experimental progress and based on the de¬ 
velopment of powerful numerical methods that allow to 
tackle these challenging problems. The nonequilibrium 
characteristics of Id systems in a homogeneous electri¬ 
cal field, resulting in a linearly decreasing potential, have 


been studied with various theoretical methods and in dif¬ 
ferent geometries and gauges. One route of research con¬ 
siders Id rings in the so-called temporal gauge, i.e. the 
ring is threaded by a magnetic flux that increases linearly 
in time. The advantage of this geometry is that one can 
use periodic boundary conditions. For non-interacting 
particles the problem can be solved analytically for arbi¬ 
trary ring size in the Keldysh formalism. For a review see 
Ref. 1221 It has been found that there is no dc-current pos¬ 
sible without additional dissipation channels. The latter 
can be introduced by attaching an infinite fermionic bath 
chain to each site of the ring. As soon as electron-electron 
interactions are included, the problem can only be solved 
by numerical means and for very small ring sizes with¬ 
out dissipative bath chains. The real-time evolution for 
a 10-site ring has been studiecP^for Hubbard interaction 
by the Cranck-Nicholson method and in Ref. |M| for the 
extended Hubbard model by the Chebyshev propagation 
method. In both cases no steady state dc current has 
been found, due to the lack of dissipation processes. In 
Ref. 123] it was shown that the current characteristics are 
always dominated by BOs. 

Another class of studies considers a finite Id Hub¬ 
bard chain (central region) with linearly decreasing on¬ 
site energies, which are attached on both ends to leads, 
represented by finite tight-binding (TB) chains of non¬ 
interacting electrons. Time-dependent density-matrix 
renormalization group (tDMRG) calculations have been 
performed for a central region of 10 sites attached to 
20 lead sites on both ends.l^S The results are in qual¬ 
itative agreement with those found for the closed ring 
geometry,!^ with the quintessence of an universal di¬ 
electric breakdown characteristic of the Mott insulating 
state. It is worth mentioning that a similar result was 
found by dynamical mean-field theory (DMFT) calcula¬ 
tions on the hyper-cubic lattice.!^ However, in all real¬ 
time evolutions the current did not reach the steady- 
state. 

Yet another possible setup to study the response of one 
dimensional lattice fermions to a homogeneous electric 
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field is an infinite Id TB-chain with linearly decreasing 
on-site energies with or without electron-electron inter¬ 
actions. In the non-interacting case the temporal gauge 
has been studied in great detail in Ref. [371 where a 
steady state was obtained by attaching fermionic dissi¬ 
pative baths to each site of the chain. The dissipation 
carries away the extra energy accumulated and is essen¬ 
tial to suppress BOs and to obtain a finite steady state 
current. Also in the two-dimensional one-band Hub¬ 
bard modeP^ an explicit coupling to fermion bath chains 
was concluded to be essential for obtaining a descrip¬ 
tion of non-equilibrium steady-states. It should be noted 
that in set-ups where a finite physical device is attached 
to infinite leads, dissipation is already provided by the 
leads.l™l 

Clearly, fermionic heat baths appear rather artificial at 
first glance. However, a closer analysis reveals that their 
impact on a physical system is qualitatively the same as 
that of phonons .1^ Recent quantum many-body studies 
on this topic showed that many physical properties in 
model systems with fermion baths behave as pred icted 
by the classical Boltzmann transport theory.l2ZE2ISl po¬ 
ther arguments in favour of fermionic baths can be found 
in Refs. 135] or 1331 where the non-equilibrium steady-state 
of photo-excited correlated electrons has been studied. 
The fermion bath model was also used earlier to describe 
the electron transport in a metallic ring,l^ where it was 
denoted as an extended Biittiker’s model. It has to be 
mentioned here, that in some systems dissipation mech¬ 
anisms do not have to be included explicitly since they 
are included anyways by means of the setup. 

In the present work, we will study an infinite Hubbard- 
chain in a homogeneous electric field. Th e dissipation will 
be described by fermionic bathd^ZiaSMI attached to each 
site of the physical chain. This approach is particularly 
convenient from the theoretical point of view. In contrast 
to previous works the non-equilibrium steady state of the 
interacting many-body system is treated by a gene raliza- 
tion of Cluster Perturbation Theory (CPT)I^^ to the 
non-equilibrium situation This allows to study the 
impact of dissipation and electron-electron interactions 
on the transport of lattice fermions in an infinite homo¬ 
geneous electric field. Furthermore, non-local effects of 
the self energy can be included via CPT. It will be shown 
later, that the inclusion of those non-local effects is cru¬ 
cial to explain the resonant structures occurring in the 
current characteristics of the correlated Hubbard chain. 

The paper is organized as follows: In Sec. |ll] we give 
an introduction to the model and its solution. The cur¬ 
rent characteristics in the non-interacting case are briefly 
summarized in Sec. [Ill A| to allow for a comprehensive 
understanding of the current characteristics expected in 
infinite Id structures. The effect of on-site Coulomb 
interactions is discussed subsequently in Sec. \niB\ To 
gain a comprehensible physical explanation of the re¬ 
sults, we describe an alternative way of modelling a non¬ 
interacting system with a similar non-equilibrium be¬ 
haviour as the interacting one in Sec. [HI C| We will close 



FIG. 1: (Color online) Sketch of the infinite one-dimensional 
correlated chain (blue squares), coupled at each site to indi¬ 
vidual Id heat bath chains (red circles). The homogeneous 
electric field S is applied parallel to the physical chain. 


our considerations by a simplified nonequilibrium Varia¬ 


tional Cluster Approach (VCA) treatment in Sec. HID 
The final conclusion is presented in Sec. |IV| 


II. THEORETICAL TREATMENT 

The setup of the model is depicted in Fig.j^ Blue 
squares represent the correlated physical sites of the ID 
Hubbard chain. Each physical site is coupled via hop¬ 
ping to fermionic bath chains. The parameters 
as depicted in the figure, stand for the hopping along 
the physical TB chain, within the bath chains, and the 
hopping between the physical system and the heat bath 
chains. The uniform electric field is applied parallel to 
the physical chain. 


A. The Hubbard-Wannier-Stark model 


The Hamiltonian corresponding to the system depicted 
in Fig. [2 is reads 

H = Ho+Hu+Hbath . (1) 

The kinetic part of the correlated chain is described by a 
TB-Hamiltonian 
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where, as usual, i^ja) denote fermionic creation (an¬ 
nihilation) operators, for site j and spin a € {tj'l-}) and 
rij = stands for the particle number operator. 
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In this paper we use t = \. In addition, we also set the 
lattice constant and the electronic charge to one. The 
on-site energies contain the linear potential due to the 
applied homogeneous electric field 8 

( 3 ) 

The electron-electron interaction present in the physi¬ 
cal chain is modelled by the Hubbard Hamiltonian 

^ “ \) “ \) 
iG 

where U represents the repulsive on-site Coulomb inter¬ 
action. 

The dissipation that is mediated by fermionic heat- 
bath chains is described by individual TB chains attached 
to each physical site. The corresponding part of the 
Hamiltonian reads 

"Hbath — — tb X] ^a+l,ja + h.C.^ + X] C ^ajcr 

OcjG 0^3<y 

( 5 ) 

3 <7 

where dljo- creation and annihilation op¬ 

erators for particles with spin a at position a within 
the bath chain, respectively, and is the corre¬ 

sponding number operator. The third term represents 
the hybridization between the site in the physical 
chain and the corresponding bath chain. As we will see 
later it is only the ratio of the two external parameters 
Ib and v that affects our results. Therefore we have used 
V = t throughout this paper. Note that the thermal bath 
chains experience the same E-field potential energy as 
the corresponding physical site. Otherwise the hybridiza¬ 
tion would violate the gauge invariance as discussed in 
App.0 The bath chains are considered to be in an equi¬ 
librium state infinitely far away from the correlated chain 
with different chemical potentials /i given by 


B. Cluster Perturbation Theory 

One way to treat properties of interacting m any-b ody 
systems is Cluster Perturbation Theory (CPT).I2^® For 
non-interacting particles, CPT always gives the exact re¬ 
sult. In the presence of electron-electron interaction CPT 
is no longer exact; it is a first order perturbation theory 
in the inter-cluster hopping parameters. The accuracy 
of CPT increases with the size of the clusters that are 
treated exactly by numerical means. 

For the model at hand, which still has a kind of trans¬ 
lational invariance, it suffices to perform an exact diago- 
nalization for one ’central cluster’. The Green’s function 
of the central cluster is then glued together iteratively by 
CPT to form the infinite system. Possible central clusters 
are depicted in Fig.[^ 
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FIG. 2: (Color online) Central clusters of size M = 2 (a), M = 
4 (b) and M = 6 (c), which in CPT are periodically strung 
together to form the infinite physical chain. The dissipation 
mechanism is represented by the wide-band result iy. 


C. Non-equilibrium CPT for the Keldysh Green’s 
functions 



In the Keldysh Green’s function formalism the current 
jmn between two adjacent lattice sites m and n for the 
aforementioned Hamiltonian can be obtained from the 
Keldysh component of the non-equilibrium Green’s func¬ 
tion: 


The CPT approach can be general ized to Keldysh 
Green’s Functions.!^ In previous studiePEH a finite de¬ 
vice with interacting electrons has been coupled to in¬ 
finite leads. Here we will generalize these ideas to infi¬ 
nite Hubbard-chains subjected to a homogeneous electric 
field. 

The single-particle Green’s functions in this general¬ 
ization are represented as four-component matrices in 
Keldysh spacd^SHU 




^ ( 
nm \ 




( 7 ) 


9(w) 


I 0 g‘^{u:)) ■ 


( 8 ) 


In the next sections we outline how the Keldysh Green’s 
function for the infinite chain for interacting electrons in 
a homogeneous electric field can be determined. 


Lower case g denotes a Green’s functions of an isolated 
cluster here, while the labels r, a and k stand for re¬ 
tarded, advanced and Keldysh component, respectively. 
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Bold letters are used for the full Keldysh Green’s func¬ 
tion given in eq. (|^. The size of the sub-matrices 5 ’', 
g°‘ and is determined by the size of the chosen clus¬ 
ter, as we will see later. Note that this notation applies 
throughout the whole paper. We will neglect the explicit 
frequency-dependence of the Green’s functions in the fol¬ 
lowing for reasons of readability and will re-introduce it 
where necessary. 

In equilibrium, the Keldysh component is related to 
the advanced/ retarded components via 


{G'') ' = ( 5 ") '-r, (13) 

which has the same form as in equilibrium. Advanced 
and retarded Green’s function are related through G“ = 
. For the Keldysh component we obtain a special 
form of the Kadanoff-Baym equatioiP^ 

= . (14) 


/ = [/-5“][2/f(a;,/r,r)-l] , (9) 

where /i?(a;,/i,T) denotes the Fermi function depending 
on the chemical potential g, and the temperature T. 


General scheme of CPT 

The basic step of the theory is the coupling of two 
initially decoupled isolated parts, which we will label by 
a,/3 e {1,2}. For this system the Green’s function is 
cluster diagonal, i.e. g°‘^ = Sai 3 g°'°‘. Now the GPT ap¬ 
proximation for the Keldysh Green’s functions G of the 
whole coupled system, written in form of a Dyson equa¬ 
tion, reads 


ga/3 


rjnCiQ. 


X ^OiOL I \ ^ ^OLCt , rriCtOc , 

Oapg + 9 *G ^ , 

/T^a Q \ 

0 • 


(lOa) 

( 10 b) 


where the matrix multiplications have to performed in 
Keldysh space (see eq. (| 8 l). T““ stands for the matrix 
that contains the parameters for the hopping processes 
that connect the two clusters a and a. In the Id case 
with only two clusters considered here, the only non-zero 
matrix elements of T““ are given by 

rT~jl,2 _ ( _ rTi2,l 

^1,M — —t — J-M,l 


where M denotes the size of the cluster. 

We are interested in the sub-matrix of G belonging to 
the first cluster only, i.e. where a = /3 = 1. Evaluating 
eq. (101 leads to 


If g belongs to a finite size cluster in equilibrium then 
the first term in curly brackets can be neglected. Insert¬ 
ing eq. 1^ into eq. (14) yields (g“) ^ — ( 5 ’’) ^ = 2i0+. 
In the applications in the present paper the second term 
will originate from semi-infinite TB-chains in the wide¬ 
band limit. It will therefore contribute a finite value for 
F^ much larger than O'*', so that the first term can be 
ignored. 


Coupling to dissipative chains 


Now we begin with the construction of the Keldysh 
Green’s function for the infinite system in an electric 
field. We start out with the computation of the Keldysh 
Green’s function g'^{uj) for the M-site central cluster 
based on Exact Diagonalization. In this approximation 
the initially uncoupled parts are in equilibrium. There¬ 
fore it suffices to compute the retarded Green’s function. 

We now use GPT to couple the central cluster to bath 
chains which are attached according to Fig.j^ We are 
only interested in Green’s functions where both in¬ 
dices a, P in eq. ( [Io| belong to the central cluster. In this 
case eq. ( 11 ) reads 


G“ = (15a) 

r := T’^bgbbrj.bc ^ 


where b denotes the bath chains, see Fig. 1. The matrix F 
describes the dissipation and de-phasing of the transport 
along the physical chain. Due to the nearest neighbour 
hopping along the bath chains the sub-matrices of 
g^^ , where C stands for retarded, advanced or Keldysh, 
are diagonal, i.e 


G“ = -I-(11a) 

pll j.l2g22j,21 

A Dyson equation in Keldysh space of the form 


= diag(gJ‘^,52'‘’,---,5M) 


(16) 


(17) 


G — g + gF G (12) 

can be solved for the advanced/ retarded and Keldysh 
component separately. For the retarded component we 
find 


Here g’’-‘’(uj, gj) is the local Green’s function of the semi¬ 
infinite bath that is connected to site j of the central 
cluster at the first contact point, ej is the on-site energy 
and gj the chemical potential in the j-th bath chain. At 
half filling we have gj = ej. The chemical potential only 
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enters the Keldysh component. For the retarded part of 
the semi-infinite TB Green’s function we hav^^ 




1 

ts 



(18) 


with 


$5 = 0 for \uj\>2tB- 

The Keldysh component follows from eq. ([^ as 
= [ g ’^-^ - [2/f(m„t) - 1]. (19) 

In the wide-band limitPE ^4 » |a;| we find 



In this case, the self-energy F entering the Dyson equa¬ 
tion (eq. @) becomes 

r = iv'^/\tb\ := . (20) 


That means that the dissipation is be described by a 
constant dissipation parameter 7. More precisely, the 
retarded part of the cluster Green’s function follows from 
that of the isolated cluster by 


for the corresponding Green’s function in which the site 
indices are restricted to the first (leftmost) cluster of the 
semi-infinite part. 

The best way to obtain an iteration equation for this 
Green’s function is given by connecting a dissipative cen¬ 
tral cluster - with its Green’s function - to the first 
cluster of the remaining semi-infinite part by GPT (see 
Fig.|3|. The coupled system then is the same as the 
original one apart form an overall energy shift ojs = £M. 
Therefore, eq. (Ill yields 


G“(a;) = G“(w) [l + 7-21 ^“(a;)] . 

( 22 ) 

This equation can be solved by considering it as an it¬ 
eration equation starting with G^^(a;) = G°'’’(a;). We ob¬ 
tain a similar equation for the left semi-infinite chain with 
ujg = —MS. Finally the two semi-infinite chains (clus¬ 
ters) can be combined again by GPT based on eq. (Ill 

Fig-Ill) ■ 
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(^cc.r) 1 ^ ^~cc.r^ ^ . ( 21 ) 


Construction of the infinite system 

Given the Green’s function of the central cluster we 
can iteratively combine these clusters to form the infi¬ 
nite chain. This can be done best by first determining a 
right/left semi-infinite ’chain’ of these clusters with the 
coupled bath chains and finally connecting them to end 
up with the desired infinite structure. 

We begin with the construction of the semi-infinite 
part on the right half of the system. Let stand 
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FIG. 4: (Color online) Decomposition of the infinite system 
into two semi-infinite parts. The zero point of the potential 
energy is chosen in the middle of the central part. 


III. RESULTS 
A. Non-interacting case 
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FIG. 3: (Color online) Iterative construction of the right semi¬ 
infinite system. The values of the potential energy ei due to 
the applied electric field are depicted below (lattice constant 
and electronic charge are set to unity). 


Now we turn to the results obtained with the Keldysh 
GPT for the model employed in the present paper. First 
we present a short summary of the key results in the 
non-interacting case and compare the description of the 
dissipative bath chains by semi-infite TB chains to the 
wide band limit approach. 

An instructive approximation is obtained for \S\,^ 

|t| and wide-band limit, where one find^^ 


3{£) 


2taq 2 ^ 



(23) 


Here q is the charge of the particles and a the lattice 
constant.lt should be re-emphasized that we use units in 
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which q = l,a = l,t = 1. This approximate result has a 
universal maximum, characterized by 


^inax — ‘2q(lV , 

t 

Jmax — • 


(24a) 

(24b) 


As emphasized before, in the non-interacting case CPT 
provides the exact result, irrespective of the cluster size. 
Therefore we here also compare our results obta ined by 
the Keldysh CPT approach to previous works 
And indeed for {7 = 0 we find that our CPT results 
for the infinite non-interacting TB-chain are in perfect 
agreement with those for the same system obtained in 
temporal gauge.^ For non-interacting particles and in 
the wide-band limit of the bath-chains it turns out that 
the coupling of the fermionic bath-chains modifies the 
cluster (advanced) self-energy by a constant damping 7 , 
as already indicated by eq. ( pT| . 

Fig. compares the current j{£) in the physical chain 
as a function of the electric field £ for the two differ¬ 
ent approaches for the dissipative bath chains. Curves 
labeled by the value for the constant damping 7 corre¬ 
spond to the wide-band limit approach, while the ones 
where dissipation is described by semi-infinite TB chains 
are labeled by the value of ts ■ The behaviour of the cur¬ 
rent is characterized by a pronounced uni-model struc¬ 
ture. While peak position and height are reasonably well 
described for all sets of parameters by the approxima¬ 
tion in eq. (24), one observes significant deviations in the 


current curve for large values of 7 and small values of 
tg respectively. However, already for 7 = 0.1 {tg = 10) 
we find an almost perfect quantitative agreement even 
up to = 1 , where the condition £" <C t is not really 
met. Since we will restrict our further investigations to 
the 7 < 0.1 regime, the wide-band limit approach chosen 
in the following is completely justified. 


eq. (23) emphasizes the importance of the heat-bath. 
Without bath chains ( 7 = 0 ) the resulting current would 
be zero for all values of as it was also observed in recent 
studies for the ICEHand 2EE11 Hubbard model. 

The linear relationship j{£) cx £ valid for weak elec¬ 
tric fields, corresponds to Ohm’s Law formulated in the 
Drude model: 


i 

£ 




-oc — 

TT F m* 


(25) 


where r denotes the lifetime of a transport electron with 
an effective mass m* and a charge e and t he TB relation¬ 
ship m* oc j and 7 cx; ^ was used.*^^*^ The expression 
essentially represents the well-known Drude DC conduc¬ 
tivity per electron. 

The initial linear region agrees with that obtained by 
linear response theory.^^ The occurrence of a maximum 
value for the current in combination with a subsequent 
negative differential conductance can be explained by the 



FIG. 5: (Color online) Current calculated according to eq. ([^ 
in the non-interacting case {U = 0). Shown with full lines 
are the results for dissipation described by a wide-band limit 
approach. In comparison the results for dissipation described 
by semi-infinite TB chains are shown with star markers. The 
approximate current gained by evaluating eq. (231 is plotted 
for comparison (black line). The position of its maximum is 
marked by a dashed line. 


increasing relevance of BOs at larger electric field inten¬ 
sities: When the electron moves to a neighbouring site 
of lower potential energy, its kinetic energy increases by 
£ (in the present units). At the same time part of the 
energy is dissipated into the heat-bath chain. The energy 
loss due to dissipation is given in terms of the effective 
damping 7 . For ^ 7 the dissipation mechanism is 
simply too weak to dissipate enough energy and conse¬ 
quently prevent the occurrence of BOs, which results in 
a continuously decreasing current. Similar results have 
been reported recently for the 2D Hubbard model.l^ 


B. Current in the Hubbard-Wannier-Stark model 

Next we want to study the impact of electronic correla¬ 
tions. For U 0 the cluster size M becomes important. 
We start the following investigations with a two-site clus¬ 
ter (M = 2). We will also present the results for larger 
cluster sizes (M = 4 and M = 6 ) before turning back to 
the two-site cluster treatment to explain the origin of the 
arising current resonances in the Hubbard-Wannier-Stark 
model. 

Fig. i shows the dependence of the current j{£) on 
the interaction strength. In all cases we used a damp¬ 
ing parameter 7 = 0.01. For U < 0.5 we find similar 
j(£)-curves as for {7 = 0, where the current maximum 
decreases slightly with increasing U. For {7 > 1 an oscil¬ 
latory behaviour sets in and a transport gap opens up. 
This gap is however much smaller then the single-particle 
gap that one obtains for £ = 0 in the density of states. 

Similar resonances have been found experimentally in 
GaAs/Al^Gai_a;As super lattices.l^ The origin of the res¬ 
onant structures in the current visible in Fig.[^ will be 
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FIG. 6 : (Color online) Current j{£) (solid line, calculated according to eq. 0 ) for different values of the interaction strength U 
= 0.5, 1.0, 2.0 and an effective damping 7 = 0.01. The U = 0 result is shown as reference in every subplot (green dashed line). 
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FIG. 7: (Color online) Current j{£) (blue solid line, calculated according to eq. 0 ) for different values of the effective damping 
7 = 0.002, 0.01, 0.05 and a fixed value for the repulsive on-site energy U = 2.0. The corresponding U = 0-result is shown as 
reference (green dashed line). 


discussed in Sec. IIII Cl 

First, however, we will study the effect of the dissi¬ 
pation strength on the current characteristics, which are 
displayed in Fig.for a fixed value of U = 2.0 and vari¬ 
ous values of 7 . There are two major effects observable in 
Fig.0 First of all, the current maximum jmax increases 
monotonically with increasing dissipation strength. This 
is in contrast to the result at t/ = 0 , where only the 
position of the maximum current is shifted according to 
eq. (24l. But it corroborates that a non-zero DC current 
is only possible when dissipation is included. The second 
effect is a smoothing of the resonances with increasing 
dissipation strength. While one finds a very spiky struc¬ 
ture in the j{E) curve for a small effective damping, the 
oscillations for larger values of 7 are rather smooth. The 
origin of this behaviour will be discussed in the following 
section. 

Now we focus on the influence of larger cluster sizes 
on the current characteristics analyzed in Fig.|^ Inde¬ 
pendent of the interaction strength the qualitative be- 
havivour does not change significantly when the number 
of sites within the CPT cluster is increased. The spiky 
structure of the current remains as well as the finite trans¬ 
port gap. Furthermore, comparing the results for M = 4 
and M = 6 it seems that in a wide region of the elec¬ 
tric field the results coincide also in a quantitative way. 


Due to the rather weak impact of the cluster size we will 
restrict our further investigations to the M = 2 case. 


C. Origin of the resonances in the 
Hubbard-Wannier-Stark model 


We will argue in this section that the origin of the res¬ 
onant structures can be traced back to the occurence of 
short-range antiferromagnetic order. Although it is well 
known that the ID Hubbar d mo del in equilibrium does 
not exhibit long-range order , 1^21121 strong short-range fluc¬ 
tuations are nevertheless present. On the short length 
scale of the cluster, where true long-range and short- 
range order can not be distinguished, this order can be 
described reasonably well by a mean-field decoupling, i.e. 
by the Stoner model. We will show in the following, that 
the positions of the resonances are in perfect agreement 
with the maxima of the eigenvalues of the WSL obtained 
in this mean-field treatment. 

The mean-field decoupled Hamiltonian is given by 

•Hi « Hf^ = j ~ “ 1) ’ (26) 

3<y 






























































FIG. 8 : (Color online) Current j{S) (blue solid line, calculated according to eq. 0 ) for different cluster sizes M = 2, 4, 6 . The 
results are shown for an interaction strength 17 = 1.0 (top row) and U = 2.0 (bottom row) with a fixed value of the effective 
damping 7 = 0.01. The corresponding U = 0 result is shown as reference (green dashed line). 


Assuming anti-ferromagnetic order leads to 

= ( 27 ) 

JO- 

where the value for the order parameter e will be ad¬ 
justed such that the Stoner model gives the same single¬ 
particle energy gap as the Hubbard model. The inset 
of Fig. 1^ shows the equilibrium Density of States (DOS) 
p{u}) for the Hubbard and the Stoner model. Obviously, 
the gap and the low-lying excitations are well described 
by the Stoner model. 

More importantly, we see that the current character¬ 
istics of the two models are qualitatively in very good 
agreement. Concerning the position of the resonances, 
we see even perfect quantitative agreement. This allows 
a transparent explanation of the resonances in the frame 
of the much simpler Stoner model. To this end we start 
out with the bare Stoner model (77'^®'^ = Hq + 77^^) 
without coupling to the bath chains. 

It is well known (see e.g. ref. I^TIl that non-interacting 
electrons in a periodic potential, which experience in ad¬ 
dition a homogeneous electric field, are trapped in local¬ 
ized states, corresponding to Bloch-oscillators. The ener¬ 
gies form Wannier-Stark ladders (WSL). In the present 
example of an alternating potential the WSL-energies 
have the form 

(28) 


where AEn denotes a possible energy correction. With¬ 
out the periodic potential (e = 0) the WSL-energies are 
simply En = n£, corresponding to AEn = 0. It is 
straightforward to determine the eigenvalues (and hence 
AEn) of 77 ''^®'" numerically for sufficiently large systems. 
The result is depicted in Fig.for e = 0.5 as function 
of ln(£). We observe that for £ ^ 1 the energy correc¬ 
tion AEn approaches the value e, the amplitude of the 
alternating potential. For £ <1 the energy correction os¬ 
cillates with extrema at as indicated in Fig.[^ Very 
surprisingly, it turns out that these field strengths corre¬ 
spond to the values of £ where the current maxima are 
observed in Fig.[^ This is illustrated in Fig.[lT]where the 
extrema of AEn compared with the position of the 
current resonances. 

In order to unravel the origin of these oscillations in 
AEn{£) it is instructive to solve the eigenvalue problem 
in first order perturbation theory. 

The eigenvectors of T^o can be expressed a^^ 

I'km) =E'7«-«^(C)c];|0) , (29) 

n 

where Ji{C) are Bessel functions of the first kind and 
C = Y' APP- we show that the first order energy 
corrections are given by 

AE)^'> = l^m) = Jo(2C) . 


En=Tl£ + AEn , 


























































9 


0.3 




0.1 


0 

0 0.2 0.4 0.6 0.8 1 

£ 

FIG. 9: (Color online) Comparison of the current characteristics j{£) obtained for the Hubbard model and its mean-field 
approximation, both for 7 = 0.01. For the Hubbard model U = 2.0 and in the mean-field approximation the corresponding 
order parameter is e = 0.3. The latter is adjusted such that both models have the same single-particle energy gap, as can be 
seen in the inset, which depicts the corresponding equilibrium DOS. 
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FIG. 10: (Color online) Energy correction Ai5„ = E„ — n£ 
for even n as function of the logarithmic field strength. A 
comparison is given between the first order correction and 
the exact value, which is obtained numerically by exact diag- 
onalization. Note that it is the same for all even sites 

and for odd n the sign is reversed. 
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FIG. 11: (Color online) Field values corresponding to the 
current resonances for the Stoner model, depicted in Fig.[^ 
compared with the positions of extrema of AEn{£) for the ex¬ 
act eigenvalues and those obtained by first order perturbation 
theory, as discussed in the text. 


The first order correction is also plotted in Fig.[T0]and 
compared with the exact correction Ai?„ := En — n£. 
Obviously, the energy correction is well described by the 
first order term and it also shows the oscillatory depen¬ 
dence on £. To understand this behaviour it is crucial 
to bear in mind that the Wannier-Stark states l^'n) are 
localized with mean position (m) = n and variance 


In perturbation theory, we can also determine the po¬ 
sition of the extrema of AEn\£) analytically 

5 d At,, At-, '■ 

i^)\s=s* ^ )|£=£* ^ ■ 


A 

= (30) 

For completeness the proof is given in App. ^ 

Consequently, for fields £ > \/2 the state is confined to 
a single site, n = 0 say, and the first order energy correc¬ 
tion is simply = e, corroborating the asymptotic 


behaviour in Fig. 10 for £ ^ 1. Then, upon lowering £ 
the wavefunction spreads out like l/£ and reaches grad¬ 
ually sites with alternating potential, which leads to an 
oscillatory behaviour of AE^^'> as observed in Fig. 10 


Let 


be the iz-th zero of Ji (x) then 


£* 


At 


(31) 


Interestingly, to first order, the positions of the max¬ 
ima are independent of e. In Fig. El these positions are 
also included. We find that the first order result is in 
reasonable agreement with the current resonances. 
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We can now invoke second order Fermi’s golden rule to 
compute the current induced by the coupling to the bath 
chains. The first order terms vanish as the coupling to 
the bath changes the number of electrons in the physi¬ 
cal chain. Brute force application of Fermi’s golden rule 
shows that the energy correction A£’„(£) play a decisive 
role and Fermi’s golden rule corroborates the above ob¬ 
servation that the current maxima occur at the extrema 
of AEr^iS). 

Now that we have shown that the short-range order is 
responsible for the current resonances, we want to an¬ 
alyze the current suppression for small £ observed in 

Fig.|^ 

For small £ the Wannier-Stark states are the wrong 
starting point. In this case we have to consider £ as 
perturbation. Then the localized states are due to the 
alternating potential barrier and have an f-independent 
localization width. In Fig.[^ a sketch of the localized 
wavefunction centred at some site ri is depicted. The £- 
field induced linear potential allows the particle to tunnel 
through the potential barrier that ends at site r 2 . The 
escape probability is proportional to the size of the filled 
area in Fig.fl^ 

For a qualitative description we approximate the left 
half of the wave function by an exponentially decreasing 
function 



FIG. 13: (Color online) Plot of Fig.[^ zoomed into the low 
electric field region. The maxima are highlighted (magenta 
star markers) and fitted according to eq. (341 (red dashed- 
dotted lines) resulting in the ht parameters Cm, depicted 
in the figures. 


The distance r 2 — ri can be approximated by the slope 
of the linear potential r 2 — ri Ri 2ej£. Hence, we have 


j{£) = . 


(34) 


tp{r) (X exp . (32) 

The quantity ^ describes the spatial extent of the local¬ 
ized wavefunction. 


Then the escape probability and in turn the current 
should be proportional to 

j{£) fx J dr . (33) 



FIG. 12: (Color online) On-site energies (solid black line) of 
the Stoner model in a homogeneous electric field with anti¬ 
ferromagnetic order parameter. The red curve sketches the 
wavefunction localized in the potential well at site ri. The 
potential well ends on the left site at position r 2 . The escape 
probability is proportional to the size of the Hlled red area. 


We consider ^ as fit parameter, which is estimated from 
the height of the resonances for small electric fields. The 
resulting envelope is depicted in Fig.[^ The envelope 
function describes the values of the current maxima al¬ 
most perfectly for small £, which confirms our descrip¬ 
tion. The size of the transport gap, defined by the field 
strength £* for which the exponential in eq. (341 is unity, 
is characterized by 


4e 

£* = - . (35) 

It is proportional to the height potential barrier and in¬ 
verse proportional to the localization length. For large £ 
the current is suppressed by the BOs in agreement with 
the non-interacting result, which is also shown in Fig. |13| 

D. Consequences/Relevance of non-local effects in 
the self-energy 

Our CPT results deviate from a recent investigation 
of the same system within DMFT.S^ This can be espe¬ 
cially seen in the linear response behaviour: while our 
calculations suggest the infinite Id chain to be in an in¬ 
sulating state, the DMFT calculations display a metallic 
linear response. This insulating state with a too large 
gap is a known drawback of CPT also known from the 
Hubbard I approximation which corresponds to using a 
single-site cluster (M = I), as can be seen in Fig.jl^ 
(red curve). This drawback can be corrected by adding 
an additional uncorrelated bath site in the sense of a 
minimal DMFT setup. The parameters of the bath site 
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can be determine^^., wi thin the Variational Cluster 
Approach (VCA)!SiHS3| pjg ^4 curve) shows the 

resulting current in the linear response region. Here, 
it is appropriate to use the parameters determined for 
the equilibrium case £ = 0. As one can see, the in¬ 
sulating behaviour has changed to a metallic one. On 
the other hand, at half filling, the Hubbard model is ex¬ 
pected to be insulati ng for any U > 0 both in one as well 
as in two dimensions .IS3EI1 This behavior is clearly missed 
when neglecting nonlocal effects in the self-energy as in 
single-site DMFT. A minimal setup to account for the 
true physical insulating solution down to C7 = 0 consists 
of a two-site cluster {M = 2) with additional bath sites, 
as can be seen from the linear response region displayed 
in Fig. These results point out the importance of 
non-local effects in the self energy, which are induced by 
short-range antiferromagnetic fluctuations. This is also 
confirmed in an early VGA workl^ according to which 
it is more effective to include physical sites in the clus¬ 
ter, thus making the self-energy more nonlocal, rather 
than bath sites, which emphasize its dynamical charac¬ 
ter. This is the reason why we have chosen in this paper 
to use large clusters rather than smaller ones with bath 
sites. Nevertheless, it would be interesting in this con¬ 
text to investigate the same infinite Id system also by 
means of a non-equilibrium cellular DMFT (CDMFT)!^ 
approach, in particular in a doped situation. 




FIG. 14: (Color online) Linear response regime/region of the 
current j{£) (calculated according to eq. 0 ) with (blue solid 
line) and without (red solid line) an additional VGA optimized 
bath site for the single-site case {M = 1) (a) and the 2-site 
case (M = 2) (b) for [7 = 2. In the single-site treatment 
the overall linear response behaviour changes from insulating 
to metallic by adding an additional bath site, while for the 
M — 2 case the insulating behaviour remains. All results are 
shown for a fixed value of the effective damping 7 = 0.01. The 
corresponding result at [7 = 0 is shown as reference (green 
dashed line). 


IV. CONCLUSIONS 

We have seen that Hubbard chains in an electric 
field couped to heat-bath chains exhibit an oscillatory/ 
resonant behaviour in the current characteristics j{S). 
The coupling to the heat-bath is essential for a non¬ 
vanishing steady-state current. The origin of the oscilla¬ 


tions has been traced back to anti-ferromagnetic fluctu¬ 
ations which we mimicked by the Stoner model. Quite 
generally we have shown that non-interacting electrons 
in an alternating potential permeated by a homogeneous 
electric field, i.e. for a Wannier-Stark model with alter¬ 
nating on-site energies, the current exhibits an oscillatory 
behaviour as function of the field strength £, which is di¬ 
rectly linked to similar oscillations in the spacing of the 
Wannier-Stark ladder. For small field strength the cur¬ 
rent is suppressed due to localization in the periodic po¬ 
tential barriers. For field strength S much larger then the 
damping parameter F, the the current is again suppressed 
due to localization or rather Bloch-oscillations, which are 
a ubiquitous feature of the Wannier-Stark model without 
dissipation. 
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Appendix A: Equivalence of the Coulomb and the 
temporal gauge in the Hubbard-Wannier-Stark 
model 

Here we will briefly outline a proof, which is closely 
related to that given in Ref. [HSl that in the Hubbard- 
Wannier-Stark model an electric field can either be 
described by a potential term, corresponding to the 
Coulomb gauge used in eq. Q, or in a time-dependent 
Peierl’s phase, corresponding to the temporal gauge.l^ 
The two approaches are equivalent in the sense that they 
differ by a unitary transformation which leaves the Li- 
ouville equation invariant. As a consequence, the steady 
state current computed via 


j oc (Al) 

will be the same. 

To prove this statement, we start out with the 
Coulomb gauge, where the electric field adds a potential 
term 


Tie =£Y^j (rij TnJ) 

3 

to the Hamiltonian T/q without A-field. The total 
Hamiltonian is Ti = 'Hq -j- TLs, where the first term also 
contains the Hubbard interaction. In the temporal gauge, 
on the other hand, there is no f-dependent potential 
term, instead the nearest neighbour hopping parameter 
is modified to 
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U+i,i —>■ ii+i,i = h+id . 

Let the Hamiltonian in temporal gauge be denoted by 
Ti. We define a time-dependent unitary operator 


U{t) = . 

It can readily be seen that Ti. = U^'HqU, where the 
Hubbard part is not modified, as the unitary operator 
only contains density operators. Let p{t) be the density 
operator defined for the system in Coulomb gauge with 
the corresponding Liouville equation 




Then the transformed density operator p = W pU ful¬ 
fils 


(j) = 

3 

= |J/(7)P-^;|J;(7)|2 = m . 



Similarly for the variance we find 


oo 

= E 


1 — — C 


(27r)2 J ^ 

= -— E - k') . 




kk' 




Integration by parts eventually yields 


which is the Liouville equation in temporal gauge. The 
initial value for t = 0 is the same for both representations, 
i.e. po = po- It should be noted that the gauge invariance 
is violated if the bath chains do not experience the same 
E-field potential as the corresponding physical site. The 
reason is the coupling term between the physical sites and 
the thermal bath. In view of the unitary transformation 
it is obvious that the current defined in eq. ( |A1[ ) is the 
same in both representations, ji+i^i = since 


• tr(pcj'_^i = 5{f;+i.i • tr(pcj'_^^_^C; ,^)} . 


a 


2 


T 

2 


2^2 

^ ■ 


Appendix C: Perturbation theory for an alternating 
potential in a Wannier-Stark model 


Here we compute the first order energy correction of 
Ho + ■ The terms of the Hamiltonian are defined in 

eq. and eq. ( [^ . 


Appendix B: Properties of the WSL wavefunction 
for a TB model 


Here we outline some of the properties of WSL states 


of eq. (291. The key element are the Bessel functions in 


the following representation 


The mean position of the WSL state I'I'm) is 


AE;W=eEl^/-m(7)|V-' 

I 

= e“™eEl'/n(7)|V-” 

n 

_ ^ \ ' i7(sin(fe)—sin(fe')) \ ' in(k—k'+Tr) 




kk' 


27r ^ ' 




= e e*™ Jo (27) . 
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